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PRELIMINARY  ANALYSIS  OF  RESULTS  OF  A  MOUNTAIN  AREA  ATMOSPHERIC 
DIFFUSION  TEST 


Atmospheric  Testing  Technical  Group* 

I.  Preface 

Following  the  continuous  advances  in  industry,  there  has 
been  an  increasing  amount  of  harmful  gases  discharged  into  the 
atmosphere  each  day.  In  order  to  guarantee  that  the  atmospheric 
environment  of  human  life  satisfies  air  quality  standards,  it  is 
necessary  on  the  one  hand  to  rationally  lay  out  newly  con- 
structed  enterprises  and  add  controls  to  the  discharge  of  sub¬ 
stances  by  already  constructed  enterprises;  on  the  other  hand, 
it  is  necessary  to  announce  in  advance  the  future  harm  levels 
of  pollutants.  All  of  these  items  urgently  demand  an  under¬ 
standing  of  the  diluting  capabilities  of  the  atmosphere  towards 
pollutants  under  different  meteorological  conditions. 

In  recent  years,  because  there  have  been  an  increasing 
number  of  factories  built  in  varied  mountain  areas  and  along 
the  coasts,  research  on  low  level  atmospheric  structures  under 
non-uniform  topographical  conditions  and  atmospheric  diffusion 
laws  has  been  getting  an  increasing  amount  of  serious  attention. 
When  comparing  mountain  areas  and  plains,  the  atmosphere  near 
stratum  structures  and  diffusion  diluting  capability  are  very 
different.  The  problems  of  how  to  fully  use  the  beneficial 
factors  of  mountain  areas,  avoid  non-advantageous  conditions 
and  reduce  the  industrial  pollution  of  factories  constructed  in 
mountain  areas  so  that  the  natural  environment,  the  atmosphere, 
is  the  responsibility  of  engineering  design  and  rational 

^Composed  of  members  of  Qinghua  University,  Central  Meteoro¬ 
logical  Department,  Atmospheric  Physics  Institute  of  Chinese 
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of  Qinghua  University  and  the  Geological  Institute  of  the 
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industrial  layout  is  a  major  task  of  atmospheric  physics 
workers . 

This  paper  used  the  results  of  diffusion  tests  of  arti¬ 
ficial  smoke  clouds  and  neutron  activated  smoke  in  a  certain 
mountain  area  during  the  spring  and  autumn  of  1975  as  well  as 
related  meteorological  observation  data,  it  also  used  differ¬ 
ent  methods  to  calculate  the  atmospheric  diffusion  parameters 
of  this  mountain  area  especially  focusing  on  the  differences 
of  the  diffusion  diluting  capabilities  of  the  pollutants  in 
the  mountain  area  and  on  the  plain  and  compared  them  with 
related  foreign  results  whereupon  useful  results  were  obtained. 

II.  General  Situation  of  the  Test  and  Sources  of  Data 

1.  Topography.  The  test  site  was  a  mountain  area  and  its 
topography  is  as  shown  in  Fig.  1.  In  the  figure,  point  A  is 
an  iron  tower  with  a  height  of  98.2  meters,  the  base  of  the 
tower  is  located  on  a  platform  372  meters  above  sea  level.  To 
its  west,  north  and  east  are  high  mountains  and  to  the  south 
it  faces  a  level  area.  The  tower  is  located  on  one  end  of  a 
north-south  running  valley. 


Fig.  1  Topography  of  the  atmospheric  test.  The  solid  lines 
are  contour  lines  (meters),  A  is  the  iron  tower,  JL 
is  the  photographing  point,  ab  is  the  three-dimensional 
photographic  datum  line,  if  is  the  architectural 
complex,  H  is  the  chimney. 

Key:  (1)  North;  (2)  Kilometers. 


2.  Meteorological  Observations  on  the  Tower 


The  meteorological  observation  instruments  were  installed 
on  the  movable  cantilever  arm  extending  from  the  tower.  The 
arm  length  is  twice  the  length  of  the  side  of  the  tower.  We 
used  the  electrical  registering  anemograph,  thermister 
thermometer  and  hot  line  anemometer  (separately  at  the  eight 
heights  of  2,  13.4,  22.4,  31.4,  44.9,  62.9,  80.9  and  97.4 
meters)  to  measure  the  average  wind  and  temperature  gradients. 
Observations  were  done  once  an  hour  and  the  average  time  was  20 
minutes.  The  horizontal  wind  direction  pulsating  data  was 
obtained  by  an  improved  EL-2  instantaneous  wind  direction 
anemograph  (installed  at  a  height  of  49.4  meters). 

3.  Artificial  Smoke  Cloud  Diffusion  Tests 

48  kilogram  smoke  containers  were  erected  at  different 
heights  on  the  iron  tower  and  acted  as  release  sources  of  the 
artificial  smoke  clouds.  The  amount  of  smoke  produced  by  this 
method  was  uniform  and  the  smoke  issuing  time  of  each  container 
was  8-12  minutes.  In  the  downwind  direction  of  the  flue  we 
used  the  two  methods  of  the  common  camera  and  the  ground  three- 
dimensional  camera  [1]  to  continuously  photograph  the  smoke 
cloud.  Each  cloud  could  be  continuously  photographed  10  to  20 
times  and  the  sampling  time  was  about  10  minutes. 

4.  Diffusion  Tests  of  Tracer  Elements 

Twenty  tests  were  carried  out  at  heights  of  80  meters  on  the 
iron  tower.  The  emission  of  indium  nitrate  in  each  test  was 
200-500  grams  and  the  emission  mode  was  burning  an  alcohol 
solution  of  indium  nitrate  in  a  high  pressure  blowlamp  to  form 
indium  oxide  particles  with  diameters  smaller  than  0.1mm  so  as 
to  simulate  the  diffusion  movements  of  the  discharged  gases  in 
the  atmosphere.  When  released,  they  separately  arranged  fan¬ 
shaped  specimen  points  on  the  seven  arcs  downwind  0.4,  0.8, 
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1.5,  3.0,  5.0,  10.0  and  45.0  kilometers  (the  relative  elimina¬ 
tion  point  tensor  angle  was  60°). 

After  completing  the  mountain  tests,  we  also  carried  out 
5  tests  on  level  fields  about  2  kilometers  south  of  the  tower. 
Seven  specimen  points  were  distributed  equidistantly  on  four 
arcs  (in  a  60°  fan-shape)  downwind  50,  150,  250  and  450  meters. 
Each  release  was  100  grams  of  indium  nitrate. 

III.  Use  of  the  Photographic  Method  to  Determine  the  Diffusion 
Parameters  of  the  Smoke  Cloud 

a 

1.  Method  for  Calculating  Diffusion  Parameters 

On  the  downwind  side  of  the  release  point,  we  calculated 
the  basic  relationship  of  vertical  diffusion  parameter  6 
from  the  smoke  cloud's  average  contour  for  the  smoke  cloud's 
locus  photograph  as  [2] 

y.  -  «,  jin O) 

In  the  formula,  z  is  the  distance  from  the  average  flue  axis 

c 

to  the  visible  boundary  (also  called  the  vertical  half  width) , 

z  is  the  maximum  vertical  half  width  and  e  is  the  bottom  of 
n 

the  natural  logarithm. 

By  using  the  graphic  method  to  find  transcendental  equation 
(1) ,  we  can  obtain  vertical  diffusion  parameter  6  based  on 

Z 

the  z  measured  on  a  certain  downward  distance.  By  using  the 
e 

completely  analogous  method,  we  can  obtain  diffusion  parameter 
0^  of  the  beam  wind  direction. 

2.  Major  Results 

Based  on  the  various  different  conditions,  we  will  now  separ¬ 
ately  discuss  the  results  obtained  through  observations  as 
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follows: 


(1)  Diffusion  of  Leeward  Slope  of  Airflow  Passed  the 
Mountain 

When  there  is  a  northwest  wind  passing  the  northern  ridge, 
the  smoke  source's  discharge  point  on  the  tower  is  located  in 
the  cavity  area  of  the  leeward  slope  and  its  diffusion  char¬ 
acteristic  point  is  the  flue's  inclined  lower  pressure  and 
strong  perturbation  motion.  See  line  1  in  Fig.  2a  for  the 
results  of  the  average  of  the  four  data  measured  in  the  tests. 
(Afterwards,  aside  from  special  explanations,  the  stability  is 
given  based  on  Pasquill's  [3]  method  of  categorization.) 

(2)  Diffusion  Under  Local  Circulation  Control 

a.  Diffusion  When  the  Wind  is  to  the  North 

When  there  is  no  systematic  airflow  passed  the  mountain, 
there  is  a  clear  night  and  the  main  body  of  the  cold  air 
accumulated  in  the  mountain  area  flows  out  from  the  eastern 
gorge.  The  wind  direction  of  the  smoke  cloud's  discharge  point 
on  the  tower  is  predominantly  northeastern  and  the  flue  is 
relatively  stable.  See  solid  line  2  in  Fig.  2a  for  the  mean 
results  of  the  five  data  measured  from  tests. 

b.  Diffusion  When  the  Wind  is  to  the  South 

During  9  tests  when  there  was  wind  towards  the  south  the  flue 
was  even  more  stable  than  when  there  was  a  north  wind  because  the 
airflow  flowed  from  the  level  area  to  the  mountain  area  and 
perturbation  motion  was  relatively  small.  See  lines  3  and  4  in 
Fig.  2a  for  the  mean  results. 

(3)  Diffusion  of  Opposite  Level  Area 
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In  order  to  compare  and  discuss  the  differences  of  the 
diffusion  of  different  positioned  leeward  slopes,  we  also 
photographed  and  measured  the  flue  discharge  from  the  chimney 
H  (see  Fig.  1)  of  the  opposite  level  area.  Its  northern 
section  was  an  architectural  complex  with  a  500  square  meter 
range  and  the  average  height  was  20  meters.  Its  southern 
sectioh  was  relatively  level  farmland  and  the  chimney  height 
was  37  meters. 

When  there  was  a  systematic  northwest  wind  and  the  chimney 
was  still  located  in  the  wake  flow  area  of  the  mountain  area's 
rough  air,  adding  on  the  influence  of  the  architectual  complex, 
the  swinging  of  the  flue  was  relatively  large.  See  line  1  in 
Fig.  2b  for  the  mean  results  of  4  measurements. 

When  the  airflow  blows  from  the  south  section's  level  area, 
the  flue  is  relatively  level  and  straight.  See  line  2  in 
Fig.  2b  for  the  mean  values  of  three  measurements. 


•-  iliSWR.  ■«»**«*  ******»****”*»  *** 

Fig.  2  Curve  of  ff  measured  by  photographic  method  changes 
with  the  diltance. 

(See  next  page  for  key) 
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Fig.  2  (continued) 

Key:  (a)  Mountain  area  conditions  are:  smoke,  cloud 
release  height  is  80  meters  and  mean  time  is  10  minutes; 
line  1  is  strong  northwestern  wind,  D  type,  average  of 
4  times;  line  2  is  northeastern  wind,  C-D  type,  average 
of  5  times;  line  3  is  the  wind  is  to  the  south,  type  B, 
average  of  5  times;  line  4  is  wind  is  to  the  south, 
type  C,  average  of  4  times;  O  indicates  the  value  of  the 
measured  point;  (b)  Level  topographical  conditions  are: 
chimney  height  is  37  meters  and  mean  time  is  10  minutes; 
line  1  is  northwestern  wind,  type  C,  average  of  4  times; 
line  2  is  wind  is  to  the  south,  type  B-C,  average  of  3 
times;  0  indicates  the  value  of  the  measured  point; 
(l)-(3)  Meters. 


IV.  Use  of  Wind  Direction  Pulsation  Data  to  Calculate  Diffusion 
Parameters 


1.  Calculation  of  Wind  Direction  Pulsation  Standard 
Difference 

From  statistical  analysis  of  smoke  cloud  crosswise  diffusion 
and  release  point  wind  direction  pulsation  data  [4] ,  we  can  ob¬ 
tain  a  method  which  uses  the  wind  direction  pulsation  data  to 
calculate  the  density  distribution  of  the  pollutants.  Its  wind 
direction  pulsation  standard  difference  is  indicated  as 
(  CL  )~  •  *n  this,  T  is  the  sampling  time,  s=  ■■  is  the  mean 
time,  -  is  the  ratio  of  the  Lagrangian  time  scale  and  the 

®  V 

Euler ian  time  scale,  T=  —  is  the  particle  movement  time,  x  is 

the  distance  of  the  tailwind  direction  particle  movement  and  u 

is  the  mean  wind  speed  of  the  tailwind  direction.  Thus, 

( fli  )  is  the  wind  direction  pulsation  standard  difference  on 

^  >  f  s 

the  release  point  with  the  sampling  time  being  %  and  the  mean 

time  being  s=  -  .  During  tests,  we  used  the  improved  EL-2 

As  u 

spontaneous  anemorumbometer .  Its  inertial  time  is  less  than  5 
seconds,  the  paper  moving  speed  of  the  wind  direction  recorder 
is  0.5mm/second  and  in  sampling  time  ,  it  reads  one  wind 
direction  pulsation  every  5  seconds.  We  obtained  a  new 
sequence  of  different  for  the  different  mean  times  recorded 


as 


Its  variance  is  expressed  as 


-wx- 


(2) 


In  the  formula,  is  the  horizontal  wind  direction  pulsation 

value  and  N  is  the  number  of  samples.  From  formula  (2)  we  can 
derive  the  curve  of  (  A.  )  _  which  changes  with  the  mean  time. 

We  can  also  calculate  the  various  values  with  different 

degrees  of  stability. 

2.  Results  of  Calculation  of  Values 


By  using  the  wind  direction  pulsation  data  on  the  tower 
(49.4  meters)  and  the  diffusion  test  (smoke  cloud  and  tracer 
element  indium)  data,  we  can  establish  the  following  relationship 


fi 

** 


O) 


When  we  determine  the  (f  of  a  certain  distance  x  from  the  dif- 

y  — - 

fusion  tests  as  well  as  the  curve  of  observed  t/jQ  which 

changes  with  s,  we  find  and  satisfy  the  sq  value  of 


Vi  -  61 

2 


r,s0so 


on  the  curve.  Then,  from 


t  —  *!*h  (4) 

we  can  find  the  j8  value.  During  the  entire  calculation,  sampliny 
time  2? *10  minutes  and  u  is  the  mean  wind  velocity  within  10 
minutes.  For  mountain  area  conditions,  is  obtained  with 
x«1000  meters  using  the  photographic  method;  for  the  level  area, 
G*  is  obtained  from  the  concentration  distribution  data  of 

y 

indium  when  x  is  within  450  meters.  The  calculation  results  of 
related  ft  values  are  listed  in  Tables  1  and  2. 
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•ZKK 

D-E 

C-D  j 

F 

J V-  fc) 

P 

2.*  2.«  2.9  1.2  1.6  1.4 

2.1  l.S 

2.0 

1  2.0 

Table  1  JS  values  of  mountain  areas. 

Key:  (1)  Stability  type;  (2)  Mean. 


(1) 

(2) 

•sacs 

D 

B-C 

¥  15 

p 

9.5  4.5  5.3  ! 

j  1.9  2.1  3.2  3.2  6.3  8.4 

5.2 

Table  2  values  of  level  area. 

Key:  (1)  Stability  type;  (2)  Mean. 

We  can  see  from  the  tables  that  in  the  opposite  level  area, 

/&  changes  in  the  2.1  to  8.9  range,  the  mean  is  5.2  which  is 
larger  than  the  $  value  of  4  calculated  by  Hay  and  Pasquill  [4] 
(distance  source  is  100  meters)  and  it  is  basically  the  same  as 
the  (distance  source  400-800  meters)  mean  value  of  5.4  (its 
changing  range  is  2.4-10.5)  which  we  obtained  in  diffusion  tests 
carried  out  at  an  open  coastal  mountain  area.  However,  the 
value  of  /6  =2.0  for  diverse  mountain  areas  is  lower  than  the 
above  mentioned  results.  This  shows  that  the  strength  of  the 
mountain  area' is  turbulent  flow  is  larger  than  that  of  the  level 
area  and  thus  ft  tends  to  be  small;  on  the  other  hand,  the  $ 
determined  by  & '  calculated  from  the  photographed  data  still  has 
a  certain  error  and  whether  or  not  these  results  reflect  the 
actual  situation  still  awaits  more  accurate  test  data  confirma¬ 
tion. 


3.  Relationship  of  ^  and  the  sampling  Time 

The  relationship  of  0^  and  the  sampling  time  can  be 
empirically  expressed  as 
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j _  ( it)' 

\ty/ 


(5) 


From  the  6  two  hour  data  observed  on  the  tower  in  the  mountain 

area,  we  separately  calculated  the  6  sets  of  ffr.  K't- .)  of 

&  i 

^T=10,  20,  30,  60,  90  and  120  minutes  and  afterwards  using  the 
linear  regression  analysis  method  very  easily  derived  the  p 
values.  The  mean  of  the  p  values  of  6  data  was  0.21  and  its 
changing  range  is  0.14-0.31.  The  foreign  plain  areas  [5] 
usually  take  p=0.2  which  is  the  same  as  the  results  in  this 
paper . 

4.  Results  of  Calculations  of  Cf 

y 

Based  on  the  obtained  <  )^.  and  values  mentioned 

previously,  we  can  very  conveniently  derive  the  (T  values  on 
different  downwind  distances.  See  Fig.  3  for  the  results  of 
different  wind  directions  and  different  stability  types  on 
mountain  and  level  area  topography. 

(8)§ 

is 


10' 


10« 

(1)  i  iIiE,*PI!K.  C-D  m  , 

(2)  2  M  #,  li  to  «•«*.(£) 

(3)  * 

(4)  «  JpiBjg[gtj#,¥J8}»shSjio*»«l7) 

(5)  * 

Fig.  3  Calculation  of  from  wind  direction  pulsation  data, 
(see  next  page  for  key) 


Pig.  3  (continued) 

Key:  (1)  Line  1,  mountain  area,  wind  to  the  south;  C-D 
type,  average  of  18  times;  (2)  Line  2,  mountain  area, 
wind  to  the  north,  D-E  type,  average  of  18  times; 

(3)  Line  3,  mountain  area,  wind  to  the  north,  E  type, 
average  of  4  times;  (4)  Line  4,  level  area,  wind  to 
the  north,  D  type  average  of  3  times;  (5)  Line  5,  level 
area,  wind  to  the  south,  C  type,  average  of  3  times  ; 

(6)  Average  time  10  minutes,  measured  height  49.4  meters; 

(7)  Measured  height  3  meters,  mean  time  10  minutes; 

(8)  Meters;  (9)  Meters. 


V.  Use  of  Profile  Data  to  Calculate  Diffusion  Parameters 


1.  Mean  Wind  Velocity  Profile  Characteristics 


Based  on  the  temperature  vertical  gradient,  we  divided  the 
measured  mean  wind  gradient  data  (145  times)  into  the  three 
different  stability  types  of  tempera  ture  inversion,  neutral 
and  decreasing  by  degrees.  See  Fig.  4  for  the  mean  wind 
velocity  profile  flow. 


(1)1  «**««»,•  *¥*&,** 

R. 

(2) *  55 

(3) >  »n*¥«3,  m 

MR. 

(4) 0  *<>#«.). 


Fig.  4  Profiles  of  mountain  area  wind  velocity. 

Key:  (1)  Line  1  represents  temperature  inversion, 
average  of  8  times,  wind  to  the  north;  (2)  Line  2 
represents  the  neutral  condition,  average  of  55  times; 
(3)  Line  3  represents  condition  of  decreasing  by 
degrees,  wind  to  the  south.  The  solid  line  is  the 
calculated  value;  (4)  0  is  the  measured  value  (mean 
time  is  20  minutes);  (5)  Meters;  (6)  Meters/second. 
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ft. 

o  **n»(*ft*« z"  #**)• ( 3 ) 

Fig.  5  Mean  wind  velocity  profile  under  different  wind 
direction  conditions. 

Key:  (1)  Line  1  represents  the  wind  to  the  south  in 

a  neutral  condition,  average  of  25  times;  (2)  Line  2 
represents  the  wind  to  the  north  in  a  neutral  con¬ 
dition,  average  of  30  times;  (3)  0  is  the  measured 
value  (mean  time  is  20  minutes) ;  (4)  Meters; 

(5)  Meters/second. 

We  can  see  from  profile  1  (average  of  8  times)  that  up  to  a 
height  of  100  meters,  the  mean  wind  velocity  which  follows  the 
height  distribution  basically  satisfies  the  logarithmic  linear 
law's  measured  value  and  calculated  value  and  the  two  are  very 
similar.  The  wind  direction  is  mainly  to  the  north. 

Profile  2  (average  of  55  times)  represents  the  mean  wind 
velocity  profile  of  the  neutral  layer  junction.  We  can  see  that 
up  to  a  height  of  100  meters,  the  wind  basically  satisfies  the 
logarithmic  law  and  the  calculated  and  measured  values  are  the 
same.  At  the  same  time,  it  is  easy  to  obtain  the  mean  roughness 
length  as  0.62  meters  from  the  point  of  crossing  of  the  straight 
line  and  longitudinal  axis. 

We  can  see  from  profile  3  (average  of  77  times)  that  when 


below  60  meters,  it  can  basically  satisfy  the  logarithmic  law 
and  then  the  upwind  velocity  decreases  very  slowly  with  the 
height.  However,  the  entire  layer  uses  the  logarithmic  linear 
law  for  approximation  and  the  calculated  and  measured  values 
are  quite  similar.  The  wind  direction  is  mainly  to  the  south. 

In  order  to  explain  the  effects  of  surface  roughness  changes 
on  the  mean  wind  velocity  profiles,  we  categorized  the  neutral 
layer  junction  data  according  to  the  wind  direction.  Curves  1 
and  2  in  Fig.  5  separately  represent  the  mean  wind  velocity 
profiles  of  the  wind  to  the  south  and  to  the  north.  We  can  see 
from  this  that  the  wind  velocity  profiles  are  all  formed  from 
two  logarithmic  profiles  with  different  slopes.  These  two 
different  slopes  reflect  the  effects  of  the  different  roughness 
topographies  on  it.  The  low  section  shows  the  effect  of  the 
topography  near  the  tower  on  the  airflow  and  the  high  section 
shows  noticeably  different  effects  of  more  distant  topography 
as  compared  with  the  roughness  near  the  tower.  Thus,  the 
roughness  lengths  of  the  two  sections  are  also  different.  There 
are  also  foreign  results  regarding  the  effects  of  topographical 
changes  on  the  wind  velocity  profile.  It  is  unanimously  con¬ 
sidered  that  there  is  an  interface  between  the  low  layer 
effected  by  the  topography  near  the  tower  and  the  high  layer 
effected  by  more  distant  topography  and  that  the  profile  on  this 
interface  shows  points  of  inflection.  The  data  analysis  in  this 
paper  shows  that  the  mean  height  of  the  points  of  inflection  of 
the  profiles  is  about  60  meters.  The  height  of  the  single 
profile  point  of  inflection  varies  with  the  difference  of  the 
wind  direction. 

2.  Calculation  of  Logarithmic  Linear  Profile  Parameters 
and  Results 

We  can  clearly  see  from  the  discussion  in  the  above  section 
that  under  non-neutral  conditions,  the  mean  wind  velocity  pro¬ 
files  can  be  approximated  by  the  logarithmic  linear  relationship. 


The  meteorological  observations  carried  out  on  the  coast  of 
Japan  by  Sensku  [6]  show  that  the  height  suitable  for  the 
logarithmic  linear  law  can  exceed  100  meters  and  can  sometimes 
reach  200  meters.  Aside  from  a  small  number  of  profiles  which 
do  not  coincide  with  this  law  (for  example,  among  90  non¬ 
neutral  profiles,  there  were  5  parabolic  alphabets  of  lines 
under  radiation  temperature  inversion  conditions) ,  this  paper 
calculates  the  logarithmic  linear  profile  parameters  with 
various  stability  types. 

The  near  ground  layer  similarity  theory  can  only  be 
applied  when  the  stability  and  uniformity  conditions  are  satis¬ 
fied,  but  strictly  speaking,  the  mountain  areas  in  these  tests 
were  not  satisfied.  Here,  we  will  use  some  relational  form¬ 
ulas  derived  by  the  similarity  theory  as  a  type  of  empirical 
approximation  but  we  will  not  carry  out  theoretical  discussions 
From  the  similarity  theory  of  the  near  ground  layer,  we  can 
express  the  average  wind  and  temperature  which  follow  the 
distribution  of  the  height  as  [7] 


1«-  + 

(6) 

K  \  «,  L/ 

T(«)  -  T(«.)  -  T.(ta  -  +  •  -) 

'  ••  L/ 

(7) 

In  the  formula,  u*  is  the  friction  speed,  Tw  is  the  temperature 
scale,  K  is  the  Karman  constant  (usually  0.4),  zq  is  the 
roughness  length,  L  is  the  Obukhov  length  scale  and  ^ is  a 
universal  constant. 

In  the  same  way,  we  assume  that  the  average  wind  and 
temperature  profiles  of  each  observation  can  be  written  as  the 
following  experimental  expressions 
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In  the  formulas,  Au,  B^,  A^  B^,  Cu  and  CT  are  the  measured 
values  of  each  wind  and  temperature  profile  derived  by  the 
binary  linear  regression  method.  Fw+her  through  a  con¬ 
trastive  analysis  of  formulas.  (6)-(9),  we  can  derive  the  follow¬ 
ing  relationships 
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Key:  (1)  Meters/second;  (2)  Meters. 


Based  on  the  above  method,  we  calculated  a  total  of  42  profile 
examples  in  order  to  coordinate  with  the  smoke  cloud  diffusion 
test  data.  There  have  already  been  many  conclusions  abroad 
related  to  the  d#  values  and  for  comparison  the  values  of  are 
listed  in  Table  3. 
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1.07 

1.S±1.15 

Table  3  Table  of  comparison  of  devalues  [6], 

Key:  (1)  Author;  (2)  This  paper. 

We  can  see  from  the  table  that  aside  from  the  values  given 
by  Monin-Obukhov  and  Sensku,  the  a*l. 8il. 15 (0 . 5-4 . 6)  value  of 
this  paper  is  the  smallest.  This  is  possibly  related  to  the 
mountain  area  topography's  perturbation  motion  being  large 
which  causes  the  wind  velocity  changes  to  become  small. 

3.  Use  of  Profile  Data  to  Calculate  Diffusion  Parameters 

The  vertical  standard  difference  &x  can  be  given  by  the 
following  empirical  relationship  [8]: 


15 


dk. 


» _ ^ 


In  the  formula,  tf  is  the  vertical  wind  velocity  pulsation 
standard  difference  and  g(x)  is  the  distance  function.  From 
the  results  given  by  Monin  [9] ,  can  be  expressed  as 


(12) 


In  the  formula,  F(£)  is  a  universal  constant  and  when  calcul¬ 
ating  the  bottom,  all  are  approached  using  the  logarithm 
linear  relationship,  £=  j-.  This  equation  is  also  derived  under 
stable  and  uniform  conditions.  Here,  we  will  only  use  this 
relational  formula  as  a  type  of  empirical  approximation.  Based 
on  this  processing,  if  there  is  a  very  large  difference  with 
the  test  results,  this  explains  that  this  type  of  approximation 
is  not  rational;  on  the  contrary,  this  type  of  approximation  can 
be  accepted  in  actual  work  and  thus  formula  (12)  can  be  written 
as 


L  >  o 

L<  0 

In  actual  calculations,  the  g(x)  in  formula  (11)  ic  empirically 
r  2-n 

given  as  x  and  r  is  given  as  (r=  — j--)  from  the  mean  wind 
velocity  based  on  the  following  formula  of  exponent  n  changing 
according  to  the  exponential  law 


Then 


L  >  0 
L<  0 


In  the  formulas,  u  takes  the  mean  of  the  entire  layer's  wind 

2**n 

velocity,  uA,  a,  L  and  r»  — are  already  calculated  and  given, 
z  takes  the  geometric  mean  height  of  the  observed  layer  (in 
the  actual  calculations,  z-  57.4x2  14  meters).  Then,  we 

can  derive  the  relationship  of  0*2  which  follows  downwind 
distance  x  from  the  formula  (15). 

In  order  to  test  this  method  of  calculating  o*  ,  we  compared 
the  calculated  values  and  simultaneously  test  measured  ^  of 
the  smoke  cloud  diffusion.  From  the  test  data  on  half  thickness 
h  of  the  smoke  cloud's  mean  contour,  by  using  the  following 
approximation  relationship 


we  can  very  conveniently  derive  Table  4  lists  the  data  of 

9  profile  observations  and  as  well  as  completely  corresponding 
photography  observations  of  the  smoke  cloud,  (7  of  which  were  the 
B-C  type  and  2  were  the  E  type,  including  4  day  test  data),  and 
the  results  of  4  representative  observations  (including  identical 
best  and  poorest  examples) . 
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Table  4  Comparison  of  calculation  of  ff  and  measurement  of 
smoke  cloud.  2 

Key:  (1)  Time;  (2)  Meters;  (3)  Distance  X  (meters) ; 
(4)  Stability  type;  (5)  Measurement  of  smoke; 

(6)  Calculated;  (7)  Measurement  of  smoke; 

(8)  Calculated;  (9)  Measurement  of  smoke; 

(10)  Calculated;  (11)  Measurement  of  smoke; 

(12)  Calculated. 


We  can  see  from  Table  4  that  by  using  the  similarity  theory 
method,  from  the  mean  wind  and  temperature  profile  data,  the 
derived  vertical  diffusion  parameters  and  actual  situation 
have  relatively  good  uniformity.  The  uniformity  is  even  better 
under  unstable  conditions.  In  order  to  further  explain  the 
level  of  uniformity  of  the  calculated  and  measured  values,  we 
did  statistics  on  the  ratio  of  the  calculated  and  measured 
values  of  <Tz  of  the  9  examples  (40-300  .meters) .  Results  showed 
that  76.1%  had  the  ratio  smaller  than  1.5,  87.2%  had  the  ratio 
smaller  than  2  and  those  with  ratios  larger  than  2  were  examples 
of  the  E  type.  This  explains  that  even  though  the  uniform  con¬ 
ditions  of  the  mountain  area  are  difficult  to  satisfy,  when  the 
similarity  theory  method  is  used  empirically,  it  cannot  cause 
relatively  large  errors.  We  applied  this  method  and  calculated 
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a  total  of  42  examples  and  divided  them  into  5  stability  types 
based  on  length,  dimension  and  size.  Each  mean  parameter  is 
listed  in  Table  5.  See  Fig.  6  for  results  of  calculations  of 


Table  5 


Table  of  calculations  of  parameters  of 


Key:  (1)  Stability  type;  (2)  Meters;  (3)  Stability 
type;  (4)  Number  of  times;  (5)  Meters/second; 


Key:  (1)  Meters;  (2)  Meters;  (3)  I  and  II,  wind  to  south; 
(4)  V  and  VI,  wind  to  the  north;  (5)  0  represents  cal¬ 
culated  point  value. 


VI.  Use  of  Concentration  Distribution  of  Indium  to  Calculate 
Diffusion  Parameters 

We  used  the  measured  indium  ground  concentration  data  to 
separately  calculate  the  horizontal  and  vertical  diffusion  para¬ 
meters.  The  CT  is  given  by 
z 

07) 

'  ic ,  '  zc,  ) 

In  the  formula,  is  the  sampling  point's  measured  concentration. 
CT  is  derived  by  the  corsswire  integral  concentration  method. 

Its  expression  is 


e.-j* 

i.i  V  *  e.u  * 


(1*) 


In  the  formula,  when  crosswire  integral  concentration  C  ,  source 

c 

strength  Q,  mean  wind  velocity  u  amd  effective  source  height  H 
are  all  known  quantities,  we  can  then  derive 


By  using  the  above  processing  method,  we  calculated  the  test 
data  of  3  level  topographical  conditions  with  the  wind  to  the 
south.  See  Fig.  7  for  results  of 
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Fig.  7  Diffusion  parameters  calculated  from  concentration 
distribution  of  indium. 

(see  next  page  for  key) 
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fig.  7  (continued) 

Key:  (1)  Meters;  (2)  Meters;  (3)  Line  1,  level  topo¬ 

graphy  CT„,  C  type,  average  of  3  times,  wind  to  the  south; 
(4)  Line  2,  level  topography  ff ,  C  type,  average  of  3 
times,  wind  to  the  south;  (5)£zis  mountain  area  CT  ,  E 
type,  average  of  6  times,  average  time  of  one  hour,  wind 
to  the  north;  (6)  0  represents  measured  value  of  level 
topography;  (7)  Average  time  of  30  minutes. 


For  the  20  tests  on  the  mountain  areas,  because  the  topo¬ 
graphical  fluctuations  within  800  meters  were  very  large,  the 
concentration  did  not  satisfy  the  normal  distribution.  Begin¬ 
ning  from  800  meters,  the  topographical  fluctuations  were 
relatively  small  and  their  concentrations  approached  satisfying 
the  normal  distribution.  We  only  calculated  9  E  type  (wind  to 
the  north)  weather  examples  of  on  800  and  1500  meter  arcs. 
See  Fig.  7  for  mean  results. 


VII.  Discussion  of  Results 


1.  Comparison  of  Mountain  Area  and  Plain  Diffusion 
Parameters 

The  diffusion  parameters  of  mountain  area  parameters  are 
usually  larger  than  those  of  the  plain  because  of  topographical 
fluctuations.  For  comparison,  we  compared  Pasquill-Gif ford' s 
[5]  representative  values  of  ^  and  dT  as  the  plain  conditions 
interpolated  from  the  curve  of  the  diffusion  parameters  which 
change  with  the  distance  and  the  results  of  the  same  stabilities 
given  by  the  four  methods  of  this  paper.  The  ratios  of  the  800 
meter  area  diffusion  parameters  are  listed  in  Table  6. 
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Table  6  Ratio  of  mountain  area  and  plain  diffusion  parameters*. 

Key:  (1)  Type  of  weather;  (2)  Ratio;  (3)  Method; 

(4)  Indium;  (5)  Profile  method;  (6)  Photographic 
method;  (7)  Pulsation  method;  (8)  South  wind;  (9)  North 
wind;  (10)*  The  ordinary  notation  of  P  indicates  the 
value  of  the  Pasquill  plain  condition,  the  lower 
notation  of  "level"  indicates  the  observation  result 
of  level  topography  and  the  lower  notation  of  "mountain" 
indicates  the  observation  results  of  mountain  areas; 

(11)  Level;  (12)  Mountain;  (13)  Mountain;  (14)  Level; 
(15)  Level;  (16)  Mountain;  (17)  Mountain:  (18)  Mountain; 
(19)  Mountain;  (20)  Mountain. 
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I 
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We  can  see  that  the  indium  made  from  the  level  area  and  the 

C y  and  CT^  calculated  from  the  wind  direction  pulsation  data 

are  basically  the  same  as  the  Pasquill  results  with  the  same 

stability  type.  When  comparing  the  same  type  of  weather  of 

the  mountain  area  and  Pasquill,  the  mountain  area  is  1.9-2. 5 

z 

times  that  of  the  plain  and  the  mountain  area  is  2.0-2. 5 
times  that  of  the  plain.  Therefore,  we  can  generally  consider 
that  the  diffusion  parameters  of  the  mountain  area  are 
2.0-2. 5  times  those  of  the  plain's  diffusion  parameters.  When 
the  wind  is  especially  strong,  the  topographical  perturbation 
is  even  stronger.  If  we  use  the  photographic  method,  when 
there  is  a  strong  northwest  wind  of  the  D  type,  mountain  area 
cr  is  8  times  that  of  the  plain.  Therefore,  from  these  test 
results  which  are  average  conditions,  we  can  approximately 
consider  that  if  we  use  the  diffusion  parameters  of  plain  con¬ 
ditions  to  estimate  the  concentrations  of  mountain  areas,  then 
it  is  4-6  times  the  actual  height.  Houind  [10]  carried  out 
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plateau  diffusion  tests  in  certain  mountain  lands  of  a  mountain 
area  and  the  results  of  D  type  weather  Ccalculated/Cmeasured= 

6.0  are  identical  to  those  of  this  paper.  In  the  same  way, 

Start  [11]  carried  out  separate  diffusion  tests  on  the  tops  of 
mountains  and  in  mountain  valleys  and  the  neutral  condition 
(equivalent  to  the  C-D  type)  Ccalculatea/Cmeasured=5  *ls°  is 
basically  the  same  as  the  results  in  this  paper. 

2.  Relationship  of  Mountain  Area's  Diffusion  Parameters 
and  Wind  Direction 

Because  of  the  non-uniformity  of  the  mountain  area  topo¬ 
graphy's  roughness  space  distribution,  great  differences  are 
created  in  the  diffusion  parameters  which  follow  the  direction 
of  the  airflow.  If  the  wind  is  to  the  south,  the  topography 
is  relatively  level  and  when  the  wind  is  to  the  north,  the 
upstream  is  the  topography  of  the  high  mountain  fluctuations. 
These  two  conditions  have  very  large  differences  towards  the 
perturbation  of  the  airflow.  By  comparing  lines  1  and  2  and  3 
and  4  of  Fig.  2a,  we  can  see  that  the  C-D  type  of  the  northeast 
wind  is  close  to  twice  as  large  as  the  C  type  &  of  the  south 
wind  and  when  there  is  a  strong  northwest  wind,  it  is  even 
larger.  If  there  is  a  strong  northwest  wind,  the  O'  of  the 
D  type  is  also  much  larger  than  the  B  type  when  there  is  a  south 
wind.  Given  the  same  stability  type,  it  can  reach  4-5  times 
larger.  Therefore,  the  wind  direction  is  a  very  important 
factor  in  selecting  the  mountain  area's  diffusion  parameter. 

3.  The  Mean  Wind  Velocity  Profile  Characteristics  of 
the  Mountain  Area's  Near  Ground  Layer 

The  mountain  area's  mean  wind  gradient  data  of  the  near 
ground  layer  less  than  100  meters  shows  that  when  there  are  non¬ 
neutral  conditions,  the  mean  wind  velocity  profile  can  be 
approximated  by  the  logarithmic  linear  law.  Its  Rvalue  is  much 
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smaller  than  that  of  the  plateau  and  the  mean  results  of  this 
paper  is  1.8+  1.15.  The  vertical  direction  perturbation  of 
this  topography  is  large  which  creates  mixed  uniformity. 

Neutral  conditions  basically  satisfy  the  logarithmic  distribu¬ 
tion  and  the  roughness  length  is  0.62  meters. 
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